Effect of eccentricity on binary neutron star searches in Advanced LIGO 



E. A. Huerta and Duncan A. Brown 

Department of Physics, Syracuse University, Syracuse, NY 13244, USA and 
Kavli Institute for Theoretical Physics, University of California, Santa Barbara, CA 93106. 

(Dated: January 10, 2013) 

Binary neutron stars (BNSs) are the primary source of gravitational waves for the Laser Interferometer 
Gravitational- wave Observatory (LIGO) and its international partners Virgo and KAGRA. Current BNS searches 
target field binaries whose orbits will have circularized by radiation reaction before their gravitational waves en- 
ter the Advanced LIGO sensitive band at 15 Hz. It has been suggested that a population of BNSs may form 
by n-body interactions near supermassive black holes or in globular clusters and that these systems may have 
non-negligible eccentricity in the Advanced LIGO band. We show that for BNS systems with total mass of 
2AMq (6 Mq), the effect of eccentricity e < 0.02 (0.05) is negligible and a circular search is effectual for 
these binaries. For eccentricities up to e = 0.4, we investigate the selection bias caused by neglecting eccentric- 
ity in BNS searches. If such high eccentricity systems exist, searches that specifically target eccentric binaries 
will be needed in Advanced LIGO and Virgo. 



I. INTRODUCTION 

Binary neutron star (BNS) mergers are expected to be rich 
sources of gravitational waves (GWs) Q and the progenitors 
of short gamma-ray bursts (SGRBs) f2], among other elec- 
tromagnetic counterparts |3|. These sources may be detected 
by advanced ground-based GW observatories w^ithin a 
horizon distance of 445 Mpc |7|. Observations of galactic 
double NS systems containing pulsars, which formed through 
binary evolution in the disk of the galaxy, suggest that five 
of these systems are expected to merge through GW radia- 
tion within a Hubble time EUll. Observed field BNSs have 
orbital periods of a few hours, merger times ~ 10^~^ Myr, 
low-kick velocities ^ 50 km , and are expected to enter the 
frequency band of ground-based GW detectors with negligible 
residual eccentricities | 8 1. For example the binary pulsar PSR 
1913+16 is expected to have an eccentricity e ~ 10~^ when 
its gravitational- wave frequency is 15Hz |8, lOjl^. Popula- 
tion synthesis models constrained by the double NS systems in 
the field predict a binary coalescence rate in Advanced LIGO 
(aLIGO) of 0.4-400 per year, with a likely value of 40 per 
year |7 |. 

In addition to the observed field binaries, observations by 
Swift have detected a population of SGRBs in or near ellipti- 
cal galaxies 112]. Some of these SGRBs have spatial offsets 
that cannot be attributed to the ejection of the progenitor com- 
pact binary from the host galaxy by the kick imparted to the 
system through the supernova explosions that create its com- 
ponents 1T31 . The necessary kick velocities for the progenitor 
compact binary to remain bound to the elliptical host galaxy 
and at the same time have its apogalacticon in the halo are 
very different to those observed in field BNSs |14|. In 1131 , 
Grindlay et al. propose that SGRBs observed in or near ellip- 
ticals may be produced by the mergers of BNSs formed by dy- 
namical capture in core-collapsed globular clusters. In these 
dense stellar environments, n-body interactions may lead to 
the formation of binaries that have sizable eccentricities at 
merger 1 15 |. Based on the globular-cluster population of el- 
liptical galaxies, and scaling from the BNS observed in the 
globular cluster Ml 5 |[T6ll . Grindlay et al. showed that around 
10 - 30% of SGRBs may be generated from BNSs in globular 



clusters, although the expected detection event rate of BNSs in 
globulars for ground-based antennas is a factor ~ 200 smaller 
than the expected rate for field BNSs iTTlfTT]!. 

Dynamical interactions in galactic nuclei may also lead to 
the formation of eccentric BNSs |18, 19 1. BNSs formed near 
galactic nuclei may enter the LIGO band before the maxi- 
mum eccentricity expected from Kozai oscillations is reached. 
After multiple periodic oscillations in the orbital elements, 
and before gravitational radiation becomes important, these 
BNSs may enter the LIGO band with high eccentricities 1 19 |. 
Nonetheless, the contribution of these type of sources to the 
BNS coalesce rates is likely to be small, namely: 10, 1 and 
0.1% of the low, realistic and high estimates in |7|. In addition 
to BNS mergers, stellar-mass black holes (BHs) are expected 
to segregate and form steep cusps around supermassive BHs in 
galactic nuclei. Emission of gravitational radiation during en- 
counters among these self-segregated BHs will drive the for- 
mation and coalescence of binary BHs that aLIGO may detect 
with a rate 1 — lO^yr"^. This optimistic estimate depends 
on the initial mass function of stars in galactic nuclei and the 
mass of the heavier BHs. These systems will be distinguish- 
able from other binaries because more than 90% of them are 
expected to have significant eccentricity (e > 0.9) when they 
enter LIGO band 1 18|. The studies quoted previously do not 
take into account natal kicks. In |20|, it is shown that if a BNS 
system receives a supernova kick which is similar in magni- 
tude to its circular speed, then the likelihood of undergoing an 
extremely fast coalescence is very high. Indeed, for BNS sys- 
tems that experience planar kicks, about 24% of them coalesce 
ten thousand times faster than a BNS system that stays in a cir- 
cular orbit immediately after the supernovae. Pulsar surveys 
may not detect these events because current search algorithms 
are not equipped to cope the strong Doppler smearing effects 
present in these systems 1201 . A pulsar can be detected if its 
lifespan is of the order of 10 Myrs and so highly eccentric 
systems which are driven to coalescence within IMyr or less 
constitute a currently 'missing' BNS population that aLIGO 
may detect. 

Motivated by the possibility of an additional population 
of eccentric binaries, we revise the study carried out in 1211 
to asses the importance of including eccentricity in BNS 
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searches in the advanced detector era. Eccentric BNSs may 
not contribute significantly to the total number of BNS events 
that may be detected, but this study quantifies the selection 
bias introduced by ignoring eccentricity in aLIGO searches. 
Based on results of radio surveys that have confirmed the exis- 
tence of pulsars with masses heavier than the canonical value 
of I.35M0 f22l, we carry out a large Monte Carlo simulation 
to address the efficacy of current search pipelines to detect 
eccentric BNS events in the mass range IMq — 3Mq, mak- 
ing use of the Zero Detuned High Power (ZDHP) noise curve 
for aLIGO 1231 . To simulate the eccentric BNS signals, we 
use the model introduced in |24|, which consists of a post- 
Newtonian (FN) model that was calibrated by comparison to 
a numerical relativity simulation of an eccentric, equal-mass 
binary black hole. This eccentric PN model combines the 
3PN conservative quasi-Keplearian orbit equations with the 
2PN evolution of the orbital elements to construct an adiabatic 
waveform that best matches the numerical relativity simula- 
tion. We use the implementation of this model, described in 
1211 . with the addition of an adaptive mesh-refinement method 
and an accelerated integrator. These improvements to the 
model were necessary to carry out a large scale Monte Carlo 
simulation to study the efficacy with which GW searches for 
circular binaries may be able to detect BNS signals with non- 
negligible eccentricity in the observational band of advanced 
detectors. In Section lU we describe the method used to con- 
struct a template bank of non-eccentric BNS waveforms. Sec- 
tion lllllpresents the results of our investigations. Finally, Sec- 
tion |IV| presents a summary of our findings and future direc- 
tions of work. 



II. BANK SIMULATION STUDIES 



where 5 = n(+/i) is the output of the detector, and n is the de- 
tector noise. The normalized overlap between any two wave- 
forms /ii, /i2 is 
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The time of coalesce tc and phase of coalescence are max- 
imized over in both the SNR and the overlap (281 . The maxi- 
mized overlap 
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gives the loss in signal-to-noise ratio incurred due to mis- 
matches in the waveforms 1291 . In this paper we compute 
the effectualness of non-eccentric templates to recover eccen- 
tric systems. To do this, we construct a 'bank' of templates 
which cover the space of circular BNS signals. The tem- 
plates are place on a hexagonal lattice using the method in- 
troduced in (30II3TI, so that the loss in signal-to-noise ratio 
between a signal and the nearest template is no more than 3%. 
The metric used to place the template grid is constructed us- 
ing TaylorF2 inspiral templates that include corrections to the 
orbital phase evolution at 2PN order. The waveform family 
we will use to create the template bank waveforms is Taylor 
T4 ISJIES, which includes corrections to the phase of the 
waveform to 3.5PN order 134H361 . We have confirmed that 
this bank is effectual for TaylorT4 waveforms before consid- 
ering the effect of eccentricity. We then compute the correla- 
tion or overlap between a simulated eccentric GW signal Qi^) 
and the best fitting non-eccentric template (/i^) by maximiz- 
ing the overlap over the templates in the bank. This is known 
as the Fitting Factor {J^J^) I29L and is defined as 



In this Section, we discuss the efficacy with which a tem- 
plate bank ^2?, '^Sl of non-eccentric waveforms can recover 
GW signals of eccentric BNS systems. The Fourier transform 
of a gravitational- wave signal h{t) is given by 



hU)= / e-2-/*/i(i)dt. (1) 

The noise- weighted inner product between two signals hi and 
h2 is given by 

^ L„ SM) 

where Sn{f) is the power spectral density of the detector 
noise. We use the form of the aLIGO ZDHP noise spec- 
trum and set the lower limit of the integral in Eq. (2) to 
/min = 15Hz. From this we can construct the matched-filter 
signal-to-noise ratio L27 J 

P=^, (3) 



TT= max 0{]f ,hl) . (6) 

6Gbank 

In the next section we present the results of this fitting factor 
investigation using the eccentric x-model waveforms as sig- 
nals. 



III. RESULTS 

The Monte Carlo simulation includes 40,000 points uni- 
formly distributed in individual component mass IM© < 
mi^2 < 3M0, and with a uniform distribution of eccentricity 
in the range e G [0, 0.4]. The template waveforms are gener- 
ated from an initial gravitational wave frequency of 15 Hz. In 
order to avoid edge effects, we generate the eccentric wave- 
forms using an initial gravitational wave frequency of 14Hz. 

The left panel of Figure [T] presents a cumulative histogram 
of the injected eccentric signals with less than the value 
associated on the x-axis. There are two intrinsic parameters 
that determine the effectualness of the circular bank to re- 
cover eccentric signals. The right panel of Figure [T] shows that 
the decreases for increasing values of eccentricity, as ex- 
pected. The right panel of Figure [T] also shows the correlation 
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FIG. 1. The left panel is a cumulative histogram which indicates the fraction of points in the mass-pace region where the bank of non-eccentric 
waveforms has a value of Fitting Factor (TT) less than the associated value on the x-axis. The right panel shows the as a function of the 
initial value of the eccentricity of the binary system. 



between the and eccentricity for systems with component 
masses mi^2 ^ [^Mq^SMq], and suggests that some systems 
can be recovered with ^ 0.95 in an eccentricity range 
e G [0.02,0.05]. 

Figure |2] shows that for BNSs the depends primarily 
on the initial value of the eccentricity and on the total mass 
of the system. It does not depend on additional parameters 
that define the dynamics of the orbit at early inspiral, i.e., the 
value of the anomaly (left panel of Figure |2]), and depends 
only weakly on the mass ratio of the binary. The right panel 
of Figure [2] shows that the tends to be larger for systems 
with larger total mass, since differences in GW phase evolu- 
tion between eccentric and circularized waveforms have less 
cycles to accumulate. In contrast, systems with low total mass 
live longer and the discrepancies in the waveforms accumulate 
over time, leading to larger phase discrepancies, and hence 
lower TTs. Based on these results, we conclude that using 
a template bank of circularized waveforms, and the place- 
ment bank algorithm described in Section |II| is be sufficient 
to recover signals with Mxotai ~ 2AMq {6Mq) emitted by 
ENS systems whose eccentricity at a fiducial GW frequency 
of 14Hz is e ^ 0.02 (0.05) with 0.95. 

We compare our findings with those reported in fSTl . Fig- 
ures 3-5 in [21 1 show the overlap between eccentric wave- 
forms generated with the eccentric time-domain x-model and 
two different circular PN approximants, namely, time-domain 
TaylorT4 2PN and the frequency domain TaylorF2 at 3.5PN. 
The differences in the overlaps between the eccentric model 
and the circular PN approximants are caused by the different 
PN orders of the models used to generate the waveforms. Fur- 
thermore, as discussed in 1371 , PN waveforms that are pertur- 
batively equivalent can have very different phase evolutions 
and asymptotic structure. These features are clearly exhibited 
in Figures 3-5 in [21 1. Nonetheless, as pointed out in 1211 . 
maximization over the template bank gives rise to a significant 
gain in overlap. Having said that, we can compare our results 



with those computed in the context of initial LIGO in Figure 6 
of 1 21 1 . Therein, it was found that initial LIGO may be able to 
efficiently detect BNS systems with Mxotai ^ 2AMq {6Mq) 
with residual eccentricity e ^ 0.05 (0.08) at 40Hz. The au- 
thors in 121] also probed the effect of eccentricity in the con- 
text of advanced LIGO assuming the PSD given in Eq. (29) 
of 1 21 1 and using a low-frequency cut off of lOHz. Given the 
demanding computational cost of such an analysis, they in- 
jected a few hundred signals to obtain a general picture of the 
efficacy with which Advanced LIGO could recover eccentric 
BNS signals. We have improved that initial analysis, and have 
carried out a large scale analysis assuming the expected sen- 
sitivity of advanced LIGO, a filter frequency of 15Hz, and an 
injection frequency of 14Hz. Taking into account the differ- 
ences in low-frequency cut offs and detector noise spectrum, 
we notice a very good agreement between our results in the 
right-hand panel of Figure [2] and those reported in Figure 8 
of 1 21 1 . We should point out, however, that we obtain more 
accurate results for low-mass BNSs because we have injected 
the eccentric signals at 14Hz. 



IV. CONCLUSIONS 

Compact binaries formed in core-collapsed globular clus- 
ters |[T3l or near supermassive black holes in galactic nu- 
clei |[T9l may have high residual eccentricities when emit- 
ting gravitational radiation in LIGO band, however there is 
significant uncertainty on the rate of these events. An inter- 
esting prospect for advanced detectors is their ability to de- 
tect the GW radiation emitted by BNS systems in wide pre- 
supernovae orbits that are driven to coalescence in very short 
times when receiving appropriately oriented retrograde kicks. 
These type of events would be observed as a supernova fol- 
lowed by a extremely fast coalescence that takes the form of 
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FIG. 2. The left panel shows the Fitting Factor {TT) as a function of the initial value of the eccentricity and the anomaly of the trajectory of 
the binary system. The right panel shows the correlation of the with the initial orbital eccentricity and the total mass of the binary system. 



a second consecutive catastrophic event. Since both events 
may take place in the same location within months or years 
apart, these events may be effectively detected as double su- 
pernovae. Even though highly eccentric BNSs may not con- 
tribute significantly to the total number of BNSs that may be 
detected by advanced ground-based GW antennas, this study 
is timely and relevant to shed light on the selection bias that 
advanced searches of BNSs may incur by neglecting eccen- 
tricity. 

We carried out a large scale Monte Carlo analysis to show 
that assuming that the expected noise curve of advanced LIGO 
is given by the Zero Detuned High Power PSD, BNS sys- 
tems with total mass ~ 2AMq (6M0) and initial eccentricity 
e ^ 0.02 (0.05) at a fiducial GW frequency of 14Hz could 
be recovered with FJ^^ 0.95 using a template bank of cir- 
cularized waveforms. In order to detect and study the rate of 
eccentric BNS in aLIGO, a search specifically targeting these 
systems will need to be constructed. The construction of such 
a template matched filter search, as well as the extension of 
this study to higher mass systems |[38||39l, is the subject of 
future work. 
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